Using Molecules to Measure Nuclear Spin-Dependent Parity Violation 
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Nuclear spin-dependent parity violation arises from weak interactions between electrons and nu- 
cleons, and from nuclear anapole moments. We outline a method to measure such effects, using a 
Stark- interference technique to determine the mixing between opposite-parity rotational/hyperfme 
levels of ground-state molecules. The technique is applicable to nuclei over a wide range of atomic 
number, in diatomic species that are theoretically tractable for interpretation. This should provide 
data on anapole moments of many nuclei, and on previously unmeasured neutral weak couplings. 
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Up to now, atomic parity violation (PV) experiments 
have primarily focused on the PV effect arising from 
the weak charge of the nucleus Qw Q) a nuclear spin- 
independent quantity that parameterizes the electroweak 
neutral coupling between electron axial- and nucleon 
vector-currents (A e V n ). Here we propose a highly sen- 
sitive and widely applicable technique to measure nu- 
clear spin- dependent (NSD) PV effects. Such effects arise 
primarily from two underlying causes. The first is the 
nuclear anapole moment, a P-odd magnetic moment in- 
duced by weak interactions within the nucleus, which 
couples to the spin of a penetrating electron 0]. Mea- 
surements of anapole moments can provide useful data 
on purely hadronic PV interactions [1, 0|- So far, only 
one nuclear anapole moment has been measured, in 133 Cs 
Q ■ The second source of NSD-P V is the electroweak neu- 
tral coupling between electron vector- and nucleon axial- 
currents (V e A n ). This can be parameterized by two con- 
stants, C2u,d, which describe the V e A n couplings to up 
and down quarks. These are suppressed in the Standard 
Model (SM), making C2 U ,d difficult to measure and at 
present perhaps the most poorly characterized parame- 
ters in the SM [fj. However, because of this suppression, 
even moderately precise measurements of C2u,d could be 
sensitive to new physics at TeV energy scales 

Our method to measure NSD-PV exploits the prop- 
erties of diatomic molecules [1, HI, [T3| to amplify the 
observable signals. Rotational/hyperfme (HF) levels of 
opposite parity can be mixed by NSD-PV interactions, 
and are inherently close in energy. Accessible laboratory 
magnetic fields can Zeeman-shift these levels to degener- 
acy, dramatically enhancing the state mixing. The ma- 
trix element (m.e.) of the NSD-PV interaction can be 
measured with a Stark-interference technique of demon- 
strated sensitivity [TTj . Use of ground-state molecules 
leads to enhanced resolution because of their long life- 
times [ll|, [l2|]- (Two recent papers also proposed mea- 
suring NSD-PV in the ground state HF levels of heavy 
alkali atoms [H|.) The technique is applicable to a wide 
class of molecules and hence to NSD-PV couplings to 
the variety of nuclei within them. We specifically con- 
sider diatomic molecules with a single valence electron in 
a 2 £ electronic state. These are the molecular equiva- 



lent of alkali atoms, with a simple, regular structure of 
rotational/HF levels. This makes it possible to reliably 
determine the properties of the electronic wavefunction 
necessary to relate the NSD-PV m.e. to the underlying 
physics. Enhancement of NSD-PV using Zeeman-shifted 
rotational/HF levels of molecules has been proposed be- 
fore Hil- We describe for the first time a realistic 
experimental approach to exploit this enhancement, and 
emphasize both its generality and its unprecedented sen- 
sitivity. 

We illustrate the basic idea of our method using a 
Hund's case (b) 2 £ diatomic molecule [H[ , containing 
one nucleus (jVi) of spin I = \ that couples to the va- 
lence electron via NSD-PV effects, and a second (A/2) 
with I' — 0. The lowest energy levels are described by 
the Hamiltonian H = B e N 2 +7N-S + M-S + c(Ln)(S-n), 
where N is the rotational angular momentum, S = | is 
the electron spin, and n is a unit vector along the inter- 
nuclear axis [l6[ (H = 1 throughout). In most cases of in- 
terest, the rotational constant B e is much larger than the 
spin-rotation (SR) constant 7 and the HF constants b and 
c. Hence N is a good quantum number, with eigenstates 
of energy E N w B e N(N + 1) and parity P = (-1) N . 

We use a magnetic field 3 = Hz to Zeeman-shift sub- 
levels of the N p = + and 1~ manifolds of states to near 
degeneracy. Since B e ^> 7, b, c, the magnetic field nec- 
essary to bridge the rotational energy E\ — Eq « 2B e 
is large enough to strongly decouple S from I and N. 
We write the molecular states in terms of the decou- 
pled basis \N,mN)\S,ms)\I,mi}. In general, the term 
c(I • n)(S • n) in H can lead to significant mixing of such 
states. However, for simplicity in the present discussion 
we assume c <c 7, b, so that the decoupled basis states 
are, to good approximation, the energy eigenstates. The 
Zeeman effect is dominated by the coupling to S, with 
approximate Hamiltonian [Io| Hz = — .g/i^S • 3, where 
g = —2 and [1b is the Bohr magneton. Opposite- 
parity levels \4>t{mN = 0, mi)) = |0, 0)||, \)\\, mi) and 
(m' N , m'j)) = |1, m! N )\j, — ^)\^ : m'j) are degenerate 
under H + H z when £ = S ~ B e /fi B . 

Pairs of these nearly-degenerate levels can be mixed 
by NSD-PV interactions, described by the relativistic 
Hamiltonian H' P = k'^" 1 



<5 3 (r) [2]. Here k' is a di 
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mensionless number parameterizing the strength of the 
NSD-PV interaction, Gp is the Fermi constant, a is the 
standard vector of Dirac matrices, and r is the displace- 
ment of the valence electron from j\f\. Within the sub- 
space of rotational/HF levels, the effect of H' P is de- 
scribed by the effective Hamiltonian Hp a = n'WpC. 
Here, Wp characterizes the size of the m.e. of H' p di- 
agonal in the 2 E wavefunction of the valence electron 
(in the molecule- fixed frame); the dimensionless opera- 
tor C = (n x S) • I/I encodes the angular momentum 
dependence of Hp s Q. We seek to determine k', by 
measuring the NSD-PV m.c.'s iW(m' N , m'j, mjsr, mi) = 
K r Wp(ip7 (m' N , m' I )\C\ip't' (mj\r,m/)) . Time-reversal in- 
variance ensures that iW is pure imaginary. 

Ingeneral, Wp cx Z 2 , where Z is the atomic number of 
A/"i The value of Wp can be explicitly calculated, 

using a semi-empirical method developed earlier by one 
of us [IH . Knowledge of &, c and 7 from standard spec- 
troscopic data [1 91 ] is sufficient to determine the s- and 
p-wave components of the valence electron wavefunction 
near A/i, and hence Wp. The approximations used in 
this method are expected to give systematic errors of 10- 
20%. This has been explicitly verified in two cases (BaF 
and YbF) by comparison of the semi-empirical results to 
sophisticated ab initio calculations of Wp [2(| ■ 

C is a pseudoscalar, with non-zero m.e.'s C between 
states with the same value of niF = mN+ms+mi. In our 
example, \tp^(0,+^)) can mix with |^7(+1, +|)), and 

|V T + (0,-i)) with both and |^-(0,+i)). 

Level crossings between the pairs of mixing states oc- 
cur at slightly different values of Ho, because of energy 
differences in the sublevels due to HF and SR terms in 
H. The magnetic fields for the various crossings differ by 
AS < 23 since (7, b, c) < B e . 

We calculate C as well as the electric dipole m.e. 
d = (ip~^\Dn z \^) , using the definitions (n Xl n y ,n z ) = 
(sin#cos</>, sin#sin^>, cos#) and \N, m^) = Y^ N (9, <p) (a 
spherical harmonic). Here D = Dn is the electric dipole 
moment in the molecular frame. To first order, d = 0, 
since [D, S] = and m! s 7^ ras- However, the HF and 
SR terms in H cause a small mixture of states with dif- 
ferent values of ms into the crossing levels. The resulting 
induced values of d can be calculated perturbatively and 
have typical size d ~ ijD, where r\ ~ (7, b, c)/B e <C 1. 

We have emphasized the essential simplicity of this sys- 
tem by outlining an analytic approach to determining all 
relevant quantities under some approximations. How- 
ever, we have also performed full numerical calculations 
of energies and m.e.'s for the system. This reproduces 
the analytic results in the simple case described here. It 
also allows inclusion of additional effects [l(| such as HF 
structure when I' 7^ 0; nuclear spin I > 1/2 and as- 
sociated electric quadrupole HF interactions; a G-tensor 
to reflect the anisotropy of the Zeeman interaction in 
the molecule; Zeeman interactions with nuclear and rota- 
tional magnetic moments; etc. The qualitative behavior 
of the system is independent of such complicating details. 

We measure iW with a Stark-interference method de- 



veloped for use in atomic Dy [ll|. A beam of molecules 
enters a region of magnetic field 23 w 23 . Here the 
molecules are excited by laser light tuned to resonance 
with a transition to a short-lived electronic state of defi- 
nite parity. Parity selection rules ensure that only one 
level of the nearly-degenerate ground state pair (say, 
is excited, and its population is rapidly depleted 
by optical pumping. Next, the molecules enter a re- 
gion of spatially varying electric field E = E(z)z, where 
E(z) = E Q sin 2itNz/L for < z < L (N is an integer) . In 
the rest frame of molecules with velocity v = vz, this 
is a time- varying field E(t — z/v) = Eosimut, where 
uj = 2ttNv/L. The Hamiltonian H± for the two-level 
system of near-degenerate states can be written as 



H± = 





-iW + dE 



iW + dE 
A 



(1) 



where A is the small 23-dependent detuning from exact 
degeneracy under H + Hz- The wavefunction is 



\m) = c+m+) + e - iAt c_(t)i^-> = ( 



- 1 c + 

c_ 



(2) 



with c+(0) = due to the optical pumping. Assuming 
W <C (dEo, A) uj, the Schrodinger equation yields 



c+(t) = -2ie~^ 
'At 

+i sin 



A A dE 2 
cos 1 — — I sin 

OJ 



W 
A" 



2 
dE 



ujt 



At the end of the electric field region, t = T = L/v and 
luT = Ntt, regardless of v. Here c+ has the final value 



„ i at / AT\ (W dE 

= 9,P. — S n \ i 1 



c+(T) = 2e — sin 



V 2 



(3) 



Next, the population of the initially depleted state is 
measured, e.g. by collecting laser-induced fluorescence 
from this state. This yields a signal S: 



S = N \c + (T)\< 



47V shr 



./AT 



W dEo 



dE 
u> 



where Ao is the number of molecules that would be de- 
tected in the absence of the optical pumping laser and 
electric field. We define the PV asymmetry A formed by 
reversal of the electric field Eq as: 



A 



S(+E )-S(-E ) W uj 



S{+E ) + S(-E a ) A dE ' 



(4) 



Formally, A diverges for A = 0. However, even in 
ideal conditions the uncertainty in W is limited by shot 
noise to SW = {WAV2S)- 1 = A/(V2Ao sin (4f)), 
which reaches a minimum value 5VF m j n = l/(2y/2N T) 
as A — > 0. In addition, under realistic conditions the 
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magnetic field CB will not be perfectly homogeneous, so 
that the values of A must be averaged over the ensemble 
of molecules. If A is described by a distibution of mean 
value Ao and r.m.s. deviation T, then for small detun- 
ings (A,r <C 2/T) the ensemble-averaged asymmetry is 

(A) = 2^^ A f^ r2 , which vanishes at A = 0. Near- 
minimum uncertainty 5W w (JWmin is then obtained 
when T and A satisfy T < A < 2/T. 

We consider molecules with rotational constant B e < 
2ir x 30 GHz, which require magnetic fields CBo < 2 T 
for level crossings. Nearly all atoms from periodic table 
groups 2,13,16, and 17, plus some from groups 1,3,4,14, 
and 15 appear in 2 S diatomic molecules satisfying this 
criterion [19| . All have laser-accessible electronic tran- 
sitions for optical pumping. A standard technique for 
producing beams of 2 S free radicals [2l|, H3| can yield 
molecular velocities d<5x 10 4 cm/s. With an electric 
field region L ~ 5 cm, SW « SW m i n can be achieved if 
r < 2tt x 3 kHz, corresponding to an r.m.s. deviation in 
the magnetic field ffB /®o < r/B e - 10" 7 . The desired 
field strength and homogeneity are available with com- 
mercial magnetic resonance imaging magnets. Molecular 
beam fluxes F ~ 10 10 /sr/s in the N — state have 
been achieved for 2 £ free radical species such as YbF 
and CaF [2l| . No depends on the solid angle of beam in- 
tercepted (fi), the detection efficiency (77), and the frac- 
tion of population in a single Zccman sublevel of the de- 
sired isotope (/). With realistic estimates £1 ~ 10~ 5 and 
77 ~ 5 x 10~ 2 , and a typical value / ~ 10~ 2 , we expect 
N = fr]QF ~ 50 mol/s and hence SW ~ 2ttx 80 Hz/V^, 
where r is the total integration time in seconds. 

From W (measured) and Wp (calculated), we deter- 
mine k' = K 2 + K,' a + k'q Q. Here, k' 2 arises from the 
V c A n term in electron-nucleus Z°-exchange; n' a from the 
nuclear anapole moment; and k'q from the coherent ef- 
fect of Qw and the magnetic HF interaction (23[. k'q is 
small compared to the other terms, and well-understood; 
we ignore it henceforth. In any measurement on a given 
nucleus, the effects of n' a and indistinguishable. 
However, they can be separated by measurements over 
a range of nuclei (24j, since k' 2 is independent of the nu- 
clear mass A while n' a oc A 2 / 3 . Thus, in heavy nuclei the 
anapole moment dominates the NSD-PV effect, while in 
light nuclei tree-level Z° exchange is primary. 

The values of n' a and k' 2 can be estimated [J] from a 
simple model of the nucleus, consisting of a single valence 
nucleon [v — N(P) for a neutron (proton)] with orbital 
angular momentum I around a uniform core: 

< « 1.0 x io-^aV 3 ^; 4 = c^^f. (5) 

Here, K = (1+1/2) (-1) /+1 / 2 ^, g v describes the strength 
of the hadronic PV interaction between v and the core, 
and [i v is the nucleon magnetic moment in nuclear mag- 
netons. In the SM, C 2 p,n are given at tree level by [l[ 

C 2P = -C 2N = A (1 - 4 sin 2 6> w ) /2 sa 0.05, (6) 



where #w is the weak mixing angle and A « 1.25. The 
values of gp^ can be related to a set of parameters 
describing low-energy hadronic PV interactions, e.g. the 
"DDH" set 25]. Based on current knowledge, op f=s 4 — 6 
and gN ~ —(0.2 — 1) are roughly expected 0, |4|; we take 
gp = 5 and gN = —1 for concreteness. Table U shows 
a sample of experimentally accessible cases. For each 
nucleus listed, we project that k' can be measured to 
~ 10% precision in an integration time of < 1 week. 

The proposed technique should also provide excellent 
control over systematic errors. The asymmetry A result- 
ing from reversal of Eq is also odd in A and in the sign of 
"Bp (It corresponds to the P-odd invariant (S — So) 
(TH-) From Ref. [ll|, it is known that the most troubling 
systematic errors (surviving all three reversals) arise from 
stray electric fields in combination with magnetic field 
gradients. The suppressed values of d mean that, even 
before reversal of A and S, macroscopic stray fields > 0.1 
mV/cm are needed to mimic the effect of iW, for all nu- 
clei considered. In addition, we find that the complex 
angular momentum dependence of the relevant operators 
makes the ratio C/d vary widely (but detcrministically) 
in magnitude and sign between nearby level crossings in 
the same molecule. This provides an additional, powerful 
test for systematics. 

A set of values for k' from many nuclei will have im- 
pact in both nuclear and particle physics. A global fit to 
all experimental data on hadronic PV at present yields 
large error bars and is internally inconsistent [H, 0] ■ Mea- 
surements in several nuclei, with both v = N and v = P, 
should be sufficient to determine the hadronic PV pa- 
rameters responsible for n' a , to moderate accuracy. The 
limiting uncertainty is likely to arise from the calcula- 
tions of nuclear structure needed to relate the anapole 
moment to the DDH parameters. At present, accuracy 
of ~30% appears reasonable, based on the spread in cal- 
culated values [3, 0, [26| . After several measurements in 
heavier nuclei, the residual effect of n' a in the lightest 
nuclei could hence be subtracted away with ~ 30% un- 
certainty, and C 2 determined with relative uncertainty 
5C 2 /C 2 w 0.3(«4/k 2 ). This is most favorable for C 2 n, 
since n' a is expected to be suppressed for odd-N nuclei by 
the small value of gN- We project 5Cw/Cw < 20% 
is possible. C 2 p,n can be written as linear combina- 
tions of C 2Uj d; e.g., C 2 n — 0.85C 2 d — 0.40C 2u including 
SU(3) and radiative corrections [27j ■ C 2 p.N have never 
been measured directly, and other linear combinations of 
C2u,d have experimental uncertainties of 70-300% [EUl]. 
A measurement of C 2 n would be complementary to a 
planned precise measurement of 2C 2u — C 2 d [1^] . 

In summary, we propose a new technique for measur- 
ing NSD-PV in a broad range of nuclei. We plan to ex- 
perimentally implement the method using 137 BaF as the 
first case. In the longer term, these ideas might be sig- 
nificantly extended. The sensitivity might be improved 
with new molecular beam sources that promise dramati- 
cally higher F at lower v [3(| . Alternatively, it might be 
possible to increase T using trapped molecular ions [3lj . 
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TABLE I: Data relevant to the proposed measurements, for a sample of nuclei in molecules where all necessary spectroscopic 
data is available. Most parameters are defined in the text. Superscript (m) indicates the value at the level-crossing with the 
maximum value of mjr. Nuclear shell- model quantum numbers are from Ref. [32j]; molecular data from Ref. [l9| . 



Such improvements, plus more spectral data for similar 
molecular species, could widen the list of accessible nu- 
clei. For light molecules and nuclei, ab initio electronic 
and nuclear structure calculations may be possible at ac- 
curacies better than those envisioned here. Ultimately, 
the method might extend to direct measurement of k' for 



1 H and 2 H. For diatomics, the low mass of H makes B e 
and hence 23 o prohibitively large. However, use of sim- 
ple triatomic species such as HSiO may circumvent this 
difficulty. 
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